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Abstract 

The theory of classically integrable nonlinear wave equations, and the Bethe 
Ansatz systems describing massive quantum field theories defined on an in- 
finite cylinder, are related by an important mathematical correspondence 
that still lacks a satisfactory physical interpretation. In this paper we shall 
describe this link for the case of the classical and quantum versions of the 
(Tzitzeica-)Bullough-Dodd model. 
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1 Introduction 



An interesting correspondence was discovered in pQ between the functional ap- 
proach to the spectral theory of ordinary differential equations, mainly developed 
by Sibuya, Voros and collaborators [21 |3], and a series of works by Bazhanov, 
Lukyanov and Zamolodchikov [H [6], where earlier results of Baxter [7] were 
extended to the conformal field theories governing the continuum limits of certain 
integrable models (IM) on a 2d lattice. 

It turns out that the initial observation of [I] was the first hint of a very general 
mathematical scheme - an 'ODE/IM correspondence' - involving wide classes 
of multi-parameter ordinary differential equations. Their generalized eigenvalue 
problems turn out to be constrained by the same Bethe Ansatz equations as arise 
in the conformal field theory limits of certain integrable vertex models, related to 
Lie algebras EJ QUI QH CE21 CE3 M ■ 

A precise relationship between elements of this new scheme and the VT- 
symmetric quantum mechanical models of Bender, Boettcher and Meisinger |15, 
[IB] was established in [17]. This has led to proofs of spectral reality in certain 
cases |18j , that were further generalised in [19] , and insights into the loss of reality 
in others [20] . Applications of the correspondence are also being found in various 
problems arising in condensed matter physics [21[ [22| [23] . 

A generalization of this correspondence to off-critical, or massive, quantum 
field theories was subsequently obtained by Gaiotto, Moore and Neitzke [24] fol- 
lowing a different chain of discoveries, relating to string and gauge theories. As 
a net outcome of a surprising series of mathematical connections, also bringing 
in the AdS/CFT correspondence, gluon scattering amplitudes in M = 4 Super 
Yang-Mills can be now studied using powerful tools from the theory of integrable 
models [25] . 

A further important step for the understanding of the relationship between 
the AdS/CFT-related results and the original formulation of the ODE/IM corre- 
spondence was made by Lukyanov and Zamolodchikov in [26]. They showed how 
to link the classical sinh-Gordon equation to the quantum massive sine(h)-Gordon 
model, generalizing the earlier works on the ODE/IM correspondence. In doing 
this, they brought into play all the formalism developed in more than fifty years 
of soliton theory and the study of exactly-solvable nonlinear partial differential 
equations (of which sinh-Gordon is an exemplar), in particular the role played by 
the association between a nonlinear classical field theory and certain sets of linear 
differential equations. 

In this paper we shall discuss another example of an off-critical generalization 
of the ODE/IM correspondence, which will be built starting from the definition 
of a particular model on the classical side. This model is often referred to as the 
Bullough-Dodd model [27], though the special properties of the corresponding 
partial differential equation were noticed as long ago as 1907, by Tzitzeica |28j. 
Like the sinh-Gordon model, it is an example of an affine Toda field theory, or 2d 
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Toda chain, as will be explained in £J2J The quantum version of the Bullough-Dodd 
model, also known as the Izergin-Korepin model [29], plays an important role 
in framework of massive 2D quantum field theories corresponding to integrable 
perturbations of the minimal series of conformal field theories [30]. Together 
with its purely-imaginary coupling version it has been directly studied using the 
exact S-matrix approach |31] , It is also related to the scaling limit of important 
statistical-mechanical systems such as the q-state Potts models |32j and the dilute 
A n models [33]. 

The analysis in this paper follows the steps taken in [II] and |26| . Some 
further details of the Bullough-Dodd model and the ODE/IM correspondence 
can be found in |34| . while a more comprehensive presentation encompassing the 
extension to other affine Toda field theories will be the subject of a forthcoming 
publication. 



2 The Bullough-Dodd model 

The one dimensional Toda chain is defined as: 

<9 t V = 2e (2T?fc+1_2T?fc) - 2e (2vk ~ 2vk ' 1 \ (2.1) 

where k enumerates and labels the fields. This definition can be simply modified 
to yield the (1 + l)-dimensional variant of (I2.ip [35J: 

d 2 v k - d 2 V k = 2e^ k+1 ~ 2 ^ - 2e ^ k - 2 v k - 1 ). (2.2) 

Varying the periodicity conditions and symmetries of the fields rj k , the system (|2,2p 
leads to a whole family of non linear equations which are solvable by means of the 
Inverse Scattering Transform (see for example |36|). This involves an associated 
linear system, which has to reproduce the nonlinear chain through a compatibility 
condition. It turns out that the correct choice for (12.21) is 1 35 1 : 



(2.3) 



XV = (d t + V-iXd -i\- 1 C 2 )^> =0 
TV = (d x + W + i\d - i\- 1 C 2 )y = 



where A is the spectral parameter, 



Vij — dtr]i5ij ; Wij — d x r]iSij ; Cuj — C 2 ji — Cjdi-ij 

a = e ^ + W) 



(2.4) 



and 5{j is the Kronecker symbol which, in a case of a closed chain consisting of iV 
elements, is defined as: 

I 1, when i = j mod N: 
8ij = I (2-5) 
0, in all other cases. 



2 



It is also possible to give a Lagrangian formulation for the system (\2.2\) : 

C = E Q^W - e^-^ 1 ) + (2.6) 
with /u E (t,x). Then it is straightforward to verify that: 



Typically the further condition Y^k^ k = is imposed on the fieldfl which results 
in the affine Toda field theory based on the affine Dynkin diagram. For 

N = 2, with r/ 1 = — rf = r], the equations then reduce to 

Vtt — Vxx + 4sinh(2?7) = (the sinh-Gordon equation). (2-7) 

For N = 3, and with r] 3 = — r/ 1 — rf, 

4-^. + 2 e (V-V)_ 2e (-v-*» a ) = o > 

^-r ? L + 2e(V + V)_ 2e (2^-2^) =0j 

which are the equations of motion for the ckp affine Toda field theory. The Z2 
symmetry of the ai affine Dynkin diagram allows a further symmetry to be 
imposed, namely 7] 1 = —rj 2 = rj. Then (|2.8p becomes 

m - Vxx + 2e 4?? - 2e~ 2r) = (the Bullough-Dodd equation). (2.9) 

Equation (|2.9p can also be recovered directly from the compatibility of 3 x 3 
matrix operators which are defined as |35| 

r] t —i\~ 1 e~ v 
T = d x + \ iXe' 71 -iA^e - " | (2.10) 

-i\~ l e 2r > i\e~v 



i\e 2 ^ 




i e -v 










1 


lx ) 





r\ x iA e 1 

X = 8 t +\ iXe-* iA^e-" | . (2.11) 
K i\- l e 2 ^ iXe-v 

If we were to introduce light-cone coordinates as 

z = x + t z = x-t, (2.12) 



"This condition has the effect of excluding a zero mode in the "mass spectrum" of the theory: 
ml = 4sin 2 (7rn/iV), n = 0, 1, . . . , N - 1 . 
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and redefine rj — >• rj/2, (|2.9p would become 

d z d 2V (z, z) + e - ^*'*) - e 2 ^'^ = . (2.13) 

The nonlinear equation ()2. 13[) is the starting point for the construction of objects 
which are in correspondence with the quantum world. However it is important 
to stress that, for this purpose, the coordinates z and z should be considered as 
living in an 'auxiliary space' which is not the same as the space on which the 
quantum field theory under discussion will be defined. In fact it will be better 
to adopt a Euclidean metric on this auxiliary space so that z and z are complex 
conjugates of each other, though this will often be relaxed and z and z regarded 
as independent complex variables. To pave the way to the correspondence, it is 
also convenient to adopt a modified version of (|2.13p . This will be referred to as 
the modified Bullough-Dodd (mBD) equation and is 

d z d- z n+e- v{z ' s) -p(z, s 3M )p(z, s 3M )e 2v{z ' s) = 0, (the mBD equation) (2.14) 

where M > and the function p(z, s 3M ) is defined to be 

p(z,s 3M ) = z 3M -s 3M , (2.15) 

and s will turn out to take the role of a scale parameter. Equations (|2.13p and 
(I2.14p are related by a change of variables and a redefinition of the field, since, 
defining 

w{z) = [ y/p(z')dz' and w(z) = [ ^p{z')dz', (2.16) 



and shifting the field as r] — > rj + | ln(p(z)p(z )), the mBD equation reduces to 
the Bullough-Dodd equation. The mBD model has in general no explicit rota- 
tional symmetry. Instead it has a discrete symmetry dictated by the form of the 
'potential' function p(z): 

z ^ e jM z , z^e~^z. ( 2 -17) 

The solutions of the mBD equation (|2.14D which will be relevant for the cur- 
rent analysis must respect this symmetry, be continuous at every finite nonzero z, 
z, and grow more slowly than exponentially at z, z — > oo. As in [26] for the (mod- 
ified) sinh-Gordon case, we will consider a one-parameter family of such solutions, 
characterized by a logarithmic behaviour near the origin. Before listing in detail 
the exact conditions which are required, polar coordinates can be introduced to 
show these properties more neatly: 

z = pe ^ i z = pe-**. (2.18) 

In most the following z and z will be treated as independent complex variables, 
so the solution r\ will be a function of the independent variables (p, 4>), according 
to (I2TTHI) . 

It is now possible to list the properties of the sought-after solutions: 
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(i) periodicity: 

2ir 

ri{pA+ m )=r]{pA); (2-19) 
or, better, the solutions r](p, (ft) are single-valued functions on a cone with 

3M ■ 



the apex angle J??, 



2tt 

: 0~<£+ — -, 0<p<oc; (2.20) 

3M OiW 

(ii) the solutions ry(p, 0) are real-valued for real p and (ft, and finite everywhere 
on the cone C^n. > except for the apex p = 0; 

3M 

(hi) large-p asymptotic: 

r)(p,(f>) = -2Mlnp + o(l) as p oo; (2.21) 

(iv) small-p asymptotic: 

r ] (p,<ft) = -2gln(p) + 0(l) for - 1< g < 1/2. (2.22) 

The factors — 2M and — 2g in (|2.2ip and (|2.22p have been chosen, for future 
convenience, to be consistent with [TT] and |26| . 

As in [26], one can start from equation (|2.22p and develop an expansion for 
z, z ~ of the form 



g H Z-z)+ m + J2lk(z 3Mk + Z 3Mk ) 
k=l 

-V0 _ s &M e 2rj 



where t]q 7^ are integration constants, all remaining terms omitted in the expan- 
sion being uniquely determined once these constants are given. These constants 
are not new parameters, but should be fixed by demanding the consistency of this 
expansion with the remaining conditions (i)-(iii) above. 

The expansion (|2.23p remains valid with z and z regarded as independent 
complex variables and, for later analysis, it will be useful to record the form of 
(I2.22p in the so-called light-cone limit z — > (with fixed z): 

r)~-gln(zz) + r ] a+'Y(z) (2.24) 

where 7(2) = YlkLl JkZ 3Mk - 

At this stage it is possible to define the linear problem associated to the mBD 
equation 

D* = 0, D* = 0, (2.25) 
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where D and D are components of an sl(3) connection and can be found by 
rearranging, in light-cone coordinates, the matrix operators (|2.10p with an appro- 
priate redefinition of the spectral parameter A and after the introduction of the 
potential p{z): 



D=d z + 



d =d- z + 



\ " 

' -\dzr\ 


\A- 1 p(z)e' 7 




(2.26) 



(2.27) 



It is possible to deal directly with this system of six equations on \l/ to get infor- 
mation about the solutions and their asymptotics, but it turns out to be simpler 
to deal with a reduction of (12.26p . (12.27P to two third-order differential equations. 
Moreover this step will be necessary for the continuation of the analysis and to 
show the connection between the original nonlinear problem and a previously- 
studied spectral problem. This reduction can be implemented defining a vector 
solution to equations fl!T26]) . (l237# l: 



Xh^d z (e^ip) 



( 



Aae 



-X'h^d^e-^) 



(2.28) 



\X2e r i/ 2 d- z (e r id- z (e- r i'4>)) / 



The equality between each row of the two parentheses in (12.280 must always hold, 
and will be useful to find large-p asymptotics of ()2.25p . Applying D and D to the 
first and the second vectors in (]2.28p respectively, it is simple to find the system 
of two third-order linear differential equations which constitutes a compatibility 
condition on if) and V> so that \l/ really is a solution of (12.25p . The two equations 
are: 



8 3 z i{j - ((d z ri) 2 + 2d 2 z ri)d z ^ 

dU - ({d-^f + 2d 2 zV )dS 



(A 3 p(z) 
(A- 3 p(z) 



d z 7]d 2 -q ■ 
d z rjd 2 r] 



0; 
0. 



(2.29) 
(2.30) 



The next step is to determine the p — > asymptotics of solutions to ()2.25|) . To 
do this it is convenient to focus attention on eq. (I2.29p . though the result would 
not change if eq. (I2.30P was taken as the starting point instead. This equation 
is entirely in z, except for the appearance of z in r](z,z) that, for now, can be 
considered as a simple parametric dependence. Finding an asymptotic solution 
ip{z) enables a solution ^ of (|2.25|) for p — > to be determined through ()2.28p . 



^This is simply recovered solving each of the systems (|2.25f) respect to two out of three 
components of the vector 
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Substituting the asymptotic form (|2.22p of 77 into (|2.29p and considering the p — > 
limit (actually z — > 0) brings equation (|2.29p into the form 

d 3 z i{j - \g(g + 2)d g if> + \g{g + 2)^ = 0, 



z 



and seeking solutions of the form ?/> = + . . . the problem is reduced to the 
solution of the following indicial equation: 

(p-l)[p(p-2)-g(g + 2)} = 0. 

This equation gives three different leading behaviours in the z — > limit: 

V'+ ~ z -9 , ~ 2, V- ~ z 9+2 . (2.31) 

These expressions give the leading behaviours and can be adjusted with appro- 
priate multiplicative constants; once inserted into (12.28|) in the usual p — > limit, 
they provide the three asymptotic solutions \P to (|2.25p . They are: 











( e (i<t>+9)9\ 


( " ' 








K ) 


\ e {-i<}>-0)g ) 




w 







~ , 3>b ~ 1 , IP_ ~ (p 0) (2.32) 



with A = e . In order to fix the values of the constants introduced above a 
particular symmetry condition has been considered. In fact, although rj(p,4>) is a 
single- valued function on the cone (|2.20p . the connection components ()2.26j) and 
(|2.27p are not. Instead, the linear problem (|2.25p is invariant with respect to the 
transformation 

2tt 2m 

O: 6^6+ , 6->6 , (2.33) 

v V 3M 3M' v ; 

involving the shift of the spectral parameter 9 (from now on 6 and A will be 
referred as spectral parameters interchangeably, keeping in mind that A = e e ). 
Imposing this symmetry, it is straightforward to find (|2.32|) . 

On the other hand it is also possible to give the large-p asymptotics, and 
for this it is best to use a semiclassical approximation. As will be explained 
more formally below, assuming that 9 is real, a generic solution of (|2.29[) grows 
exponentially at p — )■ 00, but [21 [TJ] there are special ('sub dominant') solutions 
that decay in the sector 

4vr , 4-7T 

< 4> < — — - . (2.34) 



3M + 3 3M + 3 

In order to build such a solution to (|2.25p . the first step is the specification of a 
p — > 00 asymptotic for ()2.29j) . just as done before for the solution in the region 
near p = 0, again remembering to substitute the large-p solution of ()2.14p . In 
this particular limit equation (|2.29p reduces to a sum of a third order derivative 
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and a 'potential' Q{z) which will fix the leading term of the subdominant WKB 
solution: 

Q(z) = e 3d (z 3M -s 3M ). (2.35) 
It follows from ([2.350 that the WKB-like solution has, for M > 1/2, the form: 

/ V M+1 

-M 



i;~ ClZ -M exp ^-——e» + f(z)), (2.36) 

where c\ is an arbitrary constant and the term f(z) in the exponential carries the 
dependence on z of tp. To retrieve a more exact form of (|2.36|) . it is convenient to 
seek a solution of ([2.300 of the form 

-M „„„ / z „-6 



1> ~ c 2 z~ M exp l-^—^e-" + g(z)j (2.37) 

and then by inserting both the results in ([2.280 the compatibility between the two 
solutions will fix these integration-generated functions. This yields: 

=M+i 7 M+1 

c 2 = Cl e 9 = ce d , f(z) = -- e~ e , g(z) = e 9 , (2.38) 

' JK 1 M + l av ; M + l ' v ; 

and the final semiclassical solution of ([2.290 (in polar coordinates) is 

V> W ~ ce-V M exp cosh(0 + z(M + l)0)J . (2.39) 

By inserting ([2.390 into the first of ([2.280 it is then straightforward to obtain 



^ w ~ B e^| 1 I exp I — - cosh(0 + i(M + 1)0) ) , p -»■ oo. (2.40) 



e \ / 2p M+1 



-itf>M 



M+l 



Since the functions ^— } form a basis in the space of solutions of linear 

problem ([2.250 . a linear relation of the following form must hold: 

= Q+(0, s)#+ + Q°(0, s)#b + Q-((9, s)#_, (2.41) 

where the coefficients Q^' , which are independent of the variables p and </>, are 
functions of the spectral parameter 6 and of s, as well as of the parameter g (this 
parameter has been temporarily omitted). These coefficients will coincide (up to 
an overall constant) with Q-functions of a 2d massive quantum field theory related 
to the Lie algebra A 2 , or sit (3). 
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3 The conformal limit 



Before proceeding with the advertised correspondence, it is useful to work a little 
more on equation (|2,'29p . The problem studied up to now has an important 
connection with a known spectral problem of a particular third-order differential 
equation. In order to build an appropriate analysis and to understand where the 
various steps come from, (|2.29p will be reduced to this spectral problem and a 
brief review of previously-obtained results will be given. 
Recall that equation ()2.29j) is of the form 

Sjfy - ((d z ri) 2 + 2d 2 zV ) d z i> + (A 3 p(z, s 3M ) - d z rjd 2 z r] - d 3 z n) r/> = 0, 

and, as was said above, z plays the role of a simple parameter, so it is possible to 
take the so-called light-cone limit z — > (or, if (|2.3U|) had been chosen as starting 
point then the light-cone limit would have involved z instead of z) in which rj 
assumes the form (|2.24p . After that, the limit z ~ s — > 0, 6 — > +oo can be taken, 
with the combinations 

q ji t 30 

x = e jT+i z ^ £ = s 3M e — (3.42) 

kept finite while 

x = e'~z^0, E = s 3M e~^+i ^ 0. (3.43) 

This limit is a particular scaling limit because z, s are sent to zero and 9 to 
infinity, but they are rearranged in combinations such that the new 'variables' 
x and E do not diverge or collapse. This is an important turning point of the 
theory, since it transforms the correspondence between an ODE and a massive 
QFT into a correspondence with a conformal field theory [TT] (which is massless). 
For this reason the limit (|3.43p can be called the conformal limit, with s (or, more 
precisely, s M+1 ) taking the role of the mass scale. 

Inserting (|3.42p in (|2.29p and taking the limits just discussed, equation (|2,29p 
reduces to 



3 -5(5 + 2) (^-^1 +p(:r.E) 



y(x,E,g) = 0, (3.44) 



where y has been used instead of tp to avoid confusion among the various solutions. 
Motivated by the results of [TTJ and [2] , a generalization of a theorem due to Sibuya 
for second-order equations can now be proposed: 

Conjecture 3.1 Equation (|3.44|) has a solution y(x, E, g) with the following prop- 
erties: 

(i) y is an entire function of (x,E) though, due to the branch point in the po- 
tential at x = 0, x must in general be considered to live on a suitable cover 
of the punctured complex plane; 
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(ii) y, y' = d x y and y" = d 2 x y admit, for M > 1/2, the asymptotic representa- 
tions 



AT 1 r M+l __S_ M+l M 1 T M+1 

y~x~ M e "+ 1 , y ~ -e , y"~x~ M e M +* , (3.45) 

as x —7- oo in £/ie sector 

|arg(x)l< 3MT^ ; (3 - 46) 
(mj u is uniquely fixed by the properties i) and ii). 

The behaviour of the solution (|3.45p follows from the more general formula 

y(x, E) ~ Q(x, £)~ 1/3 exp ^- ^ Q(x, t) 1/3 d^ , (3.47) 

where y solves a generic equation d x y + Q(x,E)y = 0. This is the analogue of 
a WKB approximation for a solution of a Schrodinger equation. It is possible 
to define rotated solutions, by analogy with [13], in order to construct bases of 
solutions for (|3.42j> . For general values of k, define 

y k (x,E,g) = Lo k y(co- k x,Lu- 3Mk E,g), (3.48) 

with 



Substituting in, y k solves 



8ly k - g(g + 2) (±d x - ^ ) //,, + < -^'pl.r. L')g,. - 0. (.'MO) 

and so for G Z, (|3.48p provides a potentially-new solution to (|3.42p . However 
for now it is convenient to leave k arbitrary since fractional values will also be 
needed. We also define the Stokes sectors 



2llk 

aTg{x) ~ 3MT3 



Some discussion of the dominance and subdominance of solutions is also required, 
since, for a third-order problem, this is a little more complicated than the second- 
order case. For this we largely follow [13J. First, the behaviour specified by 
conjecture 13.11 lies in the sectors S_3 U S_i U5i U53. Furthermore, in general, 

2 2 2 2 

there are three types of asymptotic solution for large |a;|. Aside the well-behaved 
leading term x~ M exp(— x M+1 /(M + 1)), there are also solutions that behave as 
x~ M exp(e ±7ri/ ' 3 x A//+1 /(M-|-l)). (This can be traced to the three third roots of —1, 
which are —1, e 7 ™/ 3 and e" 7 ™/ 3 ). According to the sector is being considered, either 
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one or two of these solutions tend to zero for large \x\. The subdominant solution 
in any given sector will be defined to be the one that tends to zero fastest in that 
sector (it might also be called maximally subdominant). Finally, y k , up to overall 
scalar factor, is characterized as the solution of (|3.42j) which is subdominant in 
the sector S k . The asymptotic (13.45j) along with the definition (13.48P imply 



(3.52) 



for \x\ — > oo with 

x e S k _3 U5 fc „i U5 fc+ i U«S fe+ s. (3.53) 

In the region formed by the two sectors S k+ i U S k+ 3, the three asymptotics yk, 
y k +i and y k +2 live together and comparing their forms, given by (|3.52p . they turn 
out to be independent of each other. A further proof of this fact can be built 
defining the Wronskian of the three solutions: 

W kl ,k 2 ,k 3 = W[y kl ,y k2 ,y k . 3 ] (3.54) 

where the 3x3 Wronskian is defined as 

// /' f"\ 

W[f,g,h]=det [g g' g" . (3.55) 
\h h! ti'J 

For /, g and h solving (|3.42p . W[f,g,h] is independent of x, and /, g and h are 
linearly independent if and only if W[f, g, h] ^ 0. In order to demonstrate the 
independence of y k , y k +\ and y k +2i it is helpful first to compute (ki, hi, k%) = 
(-1,0,1) 

W_i,o,i = — 3i\/3 (3.56) 
and then, using the general result 

W k 

i+b,k2+b,k3+b 

(E) = W klik2tk3 (u E), (3.57) 
we see that Wk k+i k+2 is always nonzero, thus confirming the independence of 

{yk, Vk+x, Vk+2}- 

The final step for this section is to show how to build new solutions to (|3.42p 
using the y k s. Supposing that k\ and k<i differ by an integer so that e~ 
e -2k 2 m _ e -2km^ ^ can ^ Q c h ec k ec i that the function 



-2kiiri 



ZkiM( x > E >9) = ykiVki - VkuVkv ( 3 - 58 ) 
which is actually a 2 x 2 Wronskian, provides a solution of 

d 3 x z klM - g(g + 2) (±d x - ±A z klM - e~ 2k7ri (x 3M - E)z klM . (3.59) 
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Equation (|3.59|) is the adjoint of (|3.50p and if k is shifted by a half-integer, then 
z kik2 becomes a solution of the original equation (|3.50[) 

^HMH ~ 9(9 + 2) (1ft, - 1) + e- 2fc -(x^ - i?). fcl+ifc2+i . 

(3.60) 

For \k\ — k%\ < 3, the regions (j3.53|) for k\ = &2 = k overlap, and it is possible to 
get an asymptotic for Zk lt k 2 from (|3.52p . In particular, for k = 1,2,3 we have 

Z-k/2,k/2(x,E,g) ~ 2i S m(7Tk/3)x- M e- 2cos ^ k/3) W xM+ \ x^oc. (3.61) 

Considering (|3.59p and (|3.6ip . for k = 1 it follows, by uniqueness of solutions, 
that 

z_i i (x, E, g) = iV3y(x, E, g). (3.62) 

2 ' 2 

This result is limited to the (k = l)-case, as for other values of the parameter is 
not possible to retrieve enough information to uniquely pin the function down. 
On the other hand (|3.62p is useful to determine the Bethe Ansatz equations, but 
this will be shown in the next section, starting from the original problem fl2.29|> . 
Now that the analyticity properties of solutions of (|3.42p have been studied, it is 
possible to return to (12.29j> and apply to that, with the appropriate caution, all 
the results that have been found in this section. 



4 Bethe Ansatz equations 

In this section, Bethe Ansatz equations will be built starting from certain func- 
tional relations involving the coefficients Q ± already introduced at the end of £}2j 
In order to get there it is necessary to apply the same analysis described in $3] to 
reproduce similar results for the original equation (|2.29p . which is repeated below 
for convenience: 

8%if> - ((d xV f + 2^r/) d x rjj + (p(x, E) - d x rjd 2 x r] - %rj)il> = 0, (4.63) 
d%if> - {(d^) 2 + 2d 2 xV )d x ^ + (p(x, E) - d x i]dh - dlrj)i> = 0, (4.64) 

where 

x = ze 9/(M+l)^ ~ = - e -e/(M+l) E = s 3M e 38M/(M+l)^ jj, _ g 3M e -30M/(M+l) 

(4.65) 

and 

ip = ip(x,x,E,E,g), $ = ip{x,x,E,E,g). (4.66) 

First of all, the rotated solutions, analogues to (I3.48p . have to be defined. In 
order to do this, the periodicity of r](p, </>) has to be exploited since it is the only 
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information, besides the asymptotic behaviours, which is known. It turns out that 
the right form is 



tPk{x, x, E, E, g) = oj k x, u~' 6MK E, u 6MK E, g), 



-3Mk 



,3Mk i 



where 



oj = exp 



2vri 



3M + 3 



(4.67) 
(4.68) 



The function (|4.67p solves 
SgVfc " {(dxV? + 2dl V )d x ^ k + {e- 2kni p(x, E) - d^rj - c^)^ = 0, (4.69) 

where, analogously to (|3.50p . a e~ 2fc?ri -term appears. In general arbitrary integer 
or half-integer values for k will be considered, and the definition of Stokes sectors 
remains unchanged (see (|3.51|) ). Following the discussion of $3j it is straightfor- 
ward to find the \x\ — > oo behaviour of the rotated solutions (|4.67[) . (keeping in 
mind that x goes to u k x) and setting 



,e/{M+i) 



(4.70) 



^~^ M+1 )x- M exp 

i>'k — ex P 

<~ W - fc ( M+1 )x M exp 



.M+l 



-LJ 



-k(M+l) 



M + l 

„M+1 
.1 w -fc(M+l) 



1 u k(M+i) 



M+l 

X M+1 

M + l 



M + l 

z.M+1 

■' MM+l) 



M + l 



(4.71) 



-k(M+l) 



sM+1 

M + l 



(the derivatives are on x only.) 

Now, the three generic solutions {V'fc, ip k +i, ipk+2} form a basis for the solutions 
of eq. (|2.29p . Again, the proof of this property comes from the evaluation of the 
Wronskian for (^1,^2,^3) = (—1,0,1) 



(4.72) 



which shows that these three solutions are independent. 

The next step is to define the analogues of the functions (|3.58j) that, here, will 
be denoted as u kltk2 so as to avoid confusion with the variable z. They read as 

u kl ,k 2 (x,x,E,E,g) = [ip kl ip'k 2 -^faViJfo^^^jflOj (4.73) 

where, again, the z-dependence has been omitted. The solutions (|4.73|) solve 



<% u ki,k2 - {{dxrfi 2 + 2d 2 r]) d x u klM 



-2irik„ 



e -■ p(x, E) + d x r)d x r\ + d x r]j u klM = 0, 

(4.74) 

and, also in this case, it is possible to evaluate (|4.73p for k\ = —k,2 = — \ 
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,7T . 



U-k/2,k/2 ~ 2isin(-A;)x M exp 



-2 cos(— fc) 



7T , / X M+1 



+ 



M+l M + l 

which, when compared with (|4.7ip for k = 1, gives the identification 

u -l/2,l/2 = «v / 3'0- 



(4.75) 



(4.76) 



To get the Bethe Ansatz equations (BAE) eq. (14. 76ft is used, exploiting the 
relation 

U-i/2,1/2 = ip-i4>'i -ipktp' i =iVSip . (4.77) 

' ' ' 2 2 2 2 

At this point the solutions ipk can be expressed in the basis of the \z\ — > solutions 
(|2,3ip with the appropriate normalization to capture the correct behaviour in the 
conformal limit 



i> k = Q+(E, E)xt + Q° k (E, E)xl + Qt(E, £) X + 



where 
and 



Qp°(E,E) = Q^{uj- k ' iM E^ k:iM E) 



±fit, ,-k3M tp , ,k3M ; 



Xk ~ x 

xl ~ ur fe ( 9+; V +2 . 



(4.78) 
(4.79) 

(4.80) 



The functions Q ±fi (E,E), introduced in (@77SJ), are the massive off-critical ana- 
logue of the CFT-related functions Q ±fi {E) [13J, [8j. Inserting the expression 
(|4,78p into (|4.77p and considering, at first, the terms proportional to x~ 9 , a func- 
tional relation containing Q + and Q° is obtained 



iV3Q+ = Ql 1/2 <5?/2^" (9+1)/2 - gQt/2Q-i/2 ujig+1)/2 + 



Qt/2Q-i/2" {g+1)/2 + gQ± 1/2 Ql 



1/2^1/2 



(4.81) 



(g + 1) {Q + _ 1/2 Q° 1/2 ^ +1)/2 - Q+ /2 Q°-i 



/2^-l/2 



(ff+l)/2 



with the appropriate to rotation factors. Equation (|4.8ip can be written in terms 
of the variable 9 = (M + 1)/(3M) In £ as: 

iV3Q + (0) = (g+1) (Q + {9 + i|)Q°(0 - if )a^ - w^Q+p - + if) 

(4.82) 

(with s and g kept constant.) 

Considering f|4.82j) evaluated at # = 9 n +m/3 and 9 = 9 n —iir/3 with (2 + (6> n ) = 
to get 



;V3Q + (0n + wr/3) = (5 + 1) (V(#„ + i27r/3)Q°(0„) 



9 + 1 
)W 2 



(4.83) 
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and 



iV3Q + (6 n - in/3) = -(g + l) L S ^Q+(9 n - i^)Q°(9 n )^J (4.84) 

respectively, and taking the ratio, the sought-after Bethe Ansatz equations are 
obtained [H] 

Q+(9 n + i2ir/3) Q+(fl n -m/3) _ {g+1) 
Q+(9 n -i2ir/3)Q+(6 n + iir/3) ' { ' ' 

The same method can be used to find BAE for the Q~ functions. In conclusion: 
by defining 9 the zeroes relative to Q^, the Bethe Ansatz equations for the two 
spectral functions are: 

Q±(9t+i2Tr/3) Q±(9$-m/Z) ±(ff+1) 

Q±(0± - z2tt/3) Q±(0± + w/3) ' l ' J 

The Bethe Ansatz equations (|4.86p can be transformed into the nonlinear integral 
equation given in section 6 of [11] for the ground state energy of the Izergin- 
Korepin model or, equivalently, of the scaling g-state Potts and tricritical Potts 
models on a cylinder geometry. This family of systems plays a key role in the 
study of integrable quantum field theories in 1+1 dimensions and statistical me- 
chanical models in 2 dimensions, and is related to the 01,2, <p2, l and 0i,5 in- 
tegrable deformations of minimal conformal field theories. The Ising model in 
external magnetic field, with its Eg-related mass spectrum [37], is perhaps the 
most physically-interesting system in the family. A schematic dictionary of the 
correspondence between the classical objects introduced here and the 'quantum 
field theory world' is given in Table [H 



classical modified Bullough-Dodd 




QFT on a cylinder 


e 




particle rapidity 


S M+1 




Rm (R: circumference, m: mass gap) 


g 




twist parameter 


M 


4$ 


l q' in the g-state Potts model 



Table 1: The dictionary 



5 Conclusions 

The main objective of this article was to discuss a particular generalization of 
the so-called ODE/IM correspondence to the massive case. The steps taken by 
Lukyanov and Zamolodchikov in |26j clarified how massive generalisations of the 
ODE/IM correspondence should be built. The results of |26j concern the relation 
between the classical and the quantum versions of the sinh-Gordon model. Adapt- 
ing their discussion it has been possible, here, to relate the classical (Titzeica-) 
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Bullough-Dodd field equation to the Izergin-Korepin massive quantum field the- 
ory. The discussion parallels precisely the conformal case treated in [11] and links 
a certain linear problem to a Bethe Ansatz system associated with the ai algebra. 
The Bethe Ansatz equations (|4.86p and the demonstration of a link with the con- 
formal field theory limit case of [H], described in section (3J are the main results 
of the paper. 

After the first discovery of this type of correspondence, it was striking how 
certain functional relations, typically emerging in the QFT domain, encoded spec- 
tral data. In the work of Lukyanov and Zamolodchikov, this rich structure was 
enlarged to contain also the relation between certain nonlinear partial differential 
equations and linear spectral problems, thereby giving further useful insights into 
the general structure of the theory. Here, further support of the validity of this 
scheme has been given, enlarging the number of working cases of the correspon- 
dence. The Bullough-Dodd equation has been chosen since, with the sinh-Gordon 
equation, it is simplest representative of the affine Toda field theories. At this 
stage it is fairly clear how a more general correspondence scheme could be devel- 
oped starting from Toda field theories based on more general Lie algebras. The 
general scheme is as follows, where the arrows have been annotated with labels of 
sections in this paper and a reference, to indicate where the corresponding step is 
further discussed for the particular case of the Bullough-Dodd equation: 



ODE 



ODE/IM correspondence 



BAE CFT 



Scaling 
limit O 



UV limit 



Linear problem 



Massive ODE/IM 
correspondence (Q 



BAE ^ massive QFT 



Zero curvature 
condition ( $2\ 



Integrable nonlinear wave equation 



Another direction for future work is the generalization of the correspondence 
between the Bethe Ansatz and classical integrable systems to non relativistically- 
invariant models such as the KdV equation and its hierarchy, and the generaliza- 
tion from integrable field theories to integrable lattice models. 

In conclusion, the connection between these two, originally disconnected, do- 
mains of mathematics and theoretical physics gives a hint of a bigger scheme in 
the wide framework of classical and quantum integrability. 
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